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Abstract. - We control the current of a single particle quantum ratchet by designing ramping 
schemes for experimentally accessible control parameters. We harvest on Landau-Zener transitions 
between Floquet states. Adiabatic and diabatic ramping allow to control the resulting directed 
transport. We find strong changes of the current in the adiabatic regime. Simple loops in control 
parameter space with alternating adiabatic and diabatic ramping are proposed. We obtain also 
current reversal. Full desymmetrization of the quantum ratchet increases the critical ramping 
speed which separates the adiabatic from the diabatic regime. 



The ratchet effect concerns the rectification of trans- 
port in the presence of zero mean forces. It was proposed 
to explain e.g. the molecular motility in biological sys- 
tems, and has been generalized to various other areas in 
physics [1,2]. Apphcations of ratchets range from biolog- 
ical systems [3] to atomic physics [4]. Especially in the 
case of light-matter interaction and laser-induced symme- 
try breaking with applications to cold atoms [4], control 
parameters are experimentally well accessible. 

A symmetry analysis of Hamiltonian and dissipative 
classical single particle ratchets was performed in Refs. 
[5,6]. A general sum rule to calculate the average veloc- 
ity of classical ratchet transport was obtained [7] . Specific 
dynamical mechanisms were also proposed [8]. Experi- 
ments on ratchet transport have also been carried out us- 
ing thermal cold rubidium and cesium atoms in optical 
lattices [4]. The time-dependent forces are applied to the 
atoms by phase modulating the laser beams which form 
the optical lattice [4] . The dissipationless case can be eas- 
ily approached by using laser beams which generate far de- 
tuned standing waves. In the latter case quantum features 
become important. To address the quantum ratchet dy- 
namics different models have been proposed mostly based 
on kicked systems [7,9, 10]. 

In two recent studies of the quantum ratchet trans- 
port [11, 12], it was found that a directed current appears 
due to a desymmetrization of the Floquet states of the sys- 



tem. A main quantum feature of the rectification process 
is the resonant interaction between Floquet states close to 
avoided crossings, which leads to a strong enhancement of 
directed transport in a narrow interval of the control pa- 
rameters. In a further extension to these studies, interac- 
tion between atoms has been considered as well [13]. In a 
different context, non-symmetric tunneling between bands 
in an asymmetric optical lattice has been reported [14]. 
The authors show that breaking the spatial symmetry in- 
duces changes of the band gap leading consequently to a 
modification of the tunneling rate between bands. 

Recent studies have also shown that it is possible to con- 
trol the states of the energy spectrum of electrons in semi- 
conductor nanostructures by ramping the external electri- 
cal field [15]. The authors exploit the diabatic and adia- 
batic Landau-Zener transitions in avoided crossings. 

We study the effect of Landau-Zener tunneling tran- 
sitions between resonant Floquet states under non- 
symmetric ac forces for a quantum ratchet. These tran- 
sitions are controlled by the ramping speed of a control 
parameter, which is related to the gap of the avoided cross- 
ing of Floquet states. We find that directed transport is 
significantly enhanced by adiabatically ramping parame- 
ters through avoided crossings between resonant Floquet 
states, even when finally ramping the systems far away 
from the resonance. We obtain current reversal, and ways 
of increasing the critical ramping velocity which separates 
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adiabatic from diabatic ramping. The results constitute 

an cxtc;nsion of the control approach in [15], from navi- 
gating between eigenstates to navigating between Floquet 
states. 

Consider a dilute gas of atoms. We can then, in good 
approximation, neglect the interaction between particles. 
The quantum dynamics for a particle moving in a periodic 
potential under the inflnonco of an ac force is described by 
the dimensionless Schrodingcr equation [11,12] 

ih^ = Hi,, (1) 

where H is the dimensionless one-particle Hamiltonian 

H=^f+ vo cos{x) - xE{t) (2) 

with p = —ih-^, and vq = 1 being the potential depth. 

ft is a dimensionless effective Planck constant. The di- 
mensionless ac field satisfies E{t + T) = E{t) with period 
T. As in [11, 12], we consider E{t) = Ei cos[w(t - to)] -f 
E2 cos[2uj{t — to) + 6] with to as initial time. By applying 
a convenient gauge transformation [11, 12, 16], we get 

H=^\p-A{t)f+cosix), (3) 

where A{t) = —Ei sin[a;(t — to)]/Lu — E2sm[2u{t — to) + 
9]/2u is the vector potential [11, 12]. 

The Hamiltonian in Eq. (2) is a periodic function of time. 
Hence the solutions of Eq.(l), \ip{t + to)) = U{to,t + 
to)\'ilj{to)), can be characterized by the eigenfunctions of 
U{to,T + to), which satisfy the relation: 

\Mt)) = e-''-^'\Mt)), \Mt + T)) = \Mt))- (4) 

The quasienergies (— tt < < tt) and the Flo- 
quet eigenstates can be obtained as solutions of the eigen- 
value problem of the Floquet operator U{O,T)\il,0{O)) = 

Duo to Bloch's theorem, discrete translational invari- 
ance of Eq.(3) implies that Floquet states arc charac- 
terized by a quasimomentum hn with \ijj(j(x + L)) = 
e^'^'^'^\ipp{x)). We choose k = (a treatment excluding 
studies of evolution of wavepackets) because our interest 
is in how a nonzero flux of atoms emerge from an initial 
state with almost zero momentum. The k — initial state 
can be produced by, for instance, the ground state of an 
elongated (shallow) trap potential of a condensate. This 
consideration allows us to use periodic boundary condi- 
tions for Eq. (1), with spatial period L = 27r. Hence the 
wave function can be expanded in the plane wave eigenba- 
sis of the momentum operator p, |n) = -^e™"^ [12]. The 
Floquet operator is obtained by solving Eqs.(l-3). Com- 
putational details are given in [12]. 

Breaking the time-reversal and shift symmetries of the 
Floquet operator [17] generates a nonzero average momen- 
tum of the Floquet states, and results in a nonzero directed 
transport [11, 12]. 



The asymptotic current of the system for t » T is ob- 
tained through the expression J{to) = X]^(l')/3|C/3(io)P 
[12], where {p)i3 are the Floquet state momenta and Cjjito) 
are the expansion coefficients of the initial wave function 
in the basis of Floquet states. The current is a func- 
tion of the relative phase 9 and the initial time to, i.e. 
J {to, 9). After averaging over the initial time the cur- 
rent becomes a function of 9 only and fulfills the relation 
J{9) = -J{9 + tt) = -J{-9) [11, 12]. 

For a given value of the control parameter 9, the ini- 
tial state \n = 0) will most strongly overlap with e.g. a 
certain Floquet state, yielding a nonzero asymptotic cur- 
rent. Variation of 9 will lead to avoided crossings between 
quasienergies, i.e. to a resonant interaction between pairs 
of Floquet states, and are confined to a small interval in 
9. Inside the avoided crossing region the Floquet states 
hybridize. If the second (new) Floquet state has a large 
mean momentum, the asymptotic current will resonantly 
increase. But upon leaving the avoided crossing region, 
the Floquet states cease to interact and to hybridize, and 
the asymptotic current returns to its nonresonant (lower) 
value. We will therefore ramp the control parameter 9 in 
time. Depending on the ramping velocity we expect to 
either reobserve the results for the case without ramping 
(diabatic Landau- Zener regime), or to slowly populate the 
second state (adiabatic regime). In the second case, af- 
ter leaving the avoided crossing region, the quantum state 
will reside mainly in the second Floquet state and keep 
the large asymptotic current. 

Hereafter the analysis is done for h = 0.2. Fig.l dis- 
plays an avoided crossing of Floquet states that leads to 
an enhancement of the current. The resonance takes place 
between a state A located in the chaotic layer of the cor- 
responding classical model, and a transporting state C 
(see [11,12] for details). We first consider an initial quan- 
tum state which exactly equals the Floquet state A in 
Fig.l. Then we ramp 9 linearly in time, viz. 

9{t) = K{9i,92,a,t) 

{9l, t < tini 

9i+ ait -tini), Uni<t<t^m + i92-9i)/a (5) 
92, t > tini + [92 - 9i)/a, 

where a is the ramping rate and tini is the time at which 
the ramping process starts. 9i and ^2 are the initial and 
final values. 

In the adiabatic regime a ^ 0, state A transforms into 
the transporting state C. The ramping time being finite, 
we estimate the deviation from the strict adiabatic behav- 
ior. We define 

i^f = a{t)i)c + b{t)i}B (6) 

as the final state in the basis of states C and B. Next we 
define [18] 

^c,B^^ dt\{^f\^|;c,B)\^ (7) 
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Fig. 1: (a) Avoided crossings between two Floquet states whidi lead to an eniiancement of tlie current (no ramping). The 
arrows indicate the direction of the ramping with initial state A. (b) Peak of the current due to the resonant states in the above 
avoided crossing, (c): Transition probability A vs a (logarithmic scale) in the basis of states B and C (see text for details). 
Circles: Numerical computation of Eq.(7). The solid lines are the fitting of Eq.(8) to the circles: P (red line), 1 — P (blue 
line).(d): Momentum evolution of the final state. The ramp is performed from 6 — —1.2 to 9 — —0.8. Momenta from top to 



bottom correspond to q = 10 , 5 x 10 ,10 
a; = 3. 
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respectively. tini = 0. The parameters are Ei = 3.26, Ei = 1.2 and 
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Fig. 2: Left panel: Current of the state |0) as a function of time in units of the period T. Ramping is performed from 6 = —1.1 
till 9 = —0.9. tini ~ 1000. Right panel: Current j vs 9 as we ramp forward and backward with different rates. We ramp 
forward from 9 = -1.1 till 9 = -0.9 with a = 10"". Then we ramp from 9 = -0.9 to 61 = -1.1 with a = 10"*, 10"*^. The 
arrows indicate the tendency of the current value as time evolves. In both panels to — 0. Parameters are the same as in Fig.l. 
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as the transition probability of the initial state A into the 

states B and C, after ramping 9 (Fig.l). 

For adiabatic ramping, we assume that variations in 
are negligible during one period of the force, which allows 
us to separate the time scales of ac driving and ramp- 
ing [18]. Then, we can use the standard Landau-Zener 
approach [19] and derive the transition probability, from 
which wc estimate the ramping rate needed for an adi- 
abatic process. The transition probability for Floquet 
states is given by [18,20], 



exp 



TT 5e 56 



2h a 



(8) 



where 66 is the 6- interval in which the quasienergy dif- 
ference increases by a factor of \/2 [18]. Se denotes the 
quasienergy gap. Thus, we find (5e « 2 x 10"'*, 66 ~ 0.01 
for the avoided crossing in Fig.l. This leads to the rough 
estimate = 5 x 10~®. To be more accurate we fit 

the transition probabilities in Fig.l using Eq.(8) and ob- 
tain « 10-'^. Therefore, for a < 10"^ the adiabatic 
regime is realized. Indeed, from Fig.lc, we see that, as the 
ramping rate decreases, the transition probability to the 
state C increases consequently. A crossover takes place 
ata~5xl0^^as the system reaches one half of the 
transition probability. 

The final state is a superposition of Floquet states B and 
C (see (6)), and therefore displays interference effects [18]. 
A signature of that interference are the oscillations expe- 
rienced by the momentum of the final state Fig. Id. These 
oscillations take place around a nonzero value, which in- 
creases as the system approaches to the adiabatic behav- 
ior. Using Eqs.(6) and (4), it is straightforward to find 

P = WPc + WPb + a*&cxp {-i^M^t) {(t)cWB) 

+ab* exp {-i^^f) {<PbWc) ■ (9) 

The last two terms in Eq.(9) generate oscillations, which 
are slow compared to the period of the driving force. The 

period of the slow oscillations in Fig. Id is in agreement 

with the estimation -. ^ r ^ t « 10^. In the extreme 

|ec - (b\ 

adiabatic limit & = 0, and consequently the last two terms 
in Eq.(9) disappear. 

Let us now address the effect of adiabatic ramping on 
the generation of directed transport. To this end, we ramp 
6 within the interval 6 G (—1.1,-0.9), where no other 
avoided crossings appear. We consider the initial state 
|0) with zero momentum. It overlaps strongly with the 
state A depicted by a circle in Fig. la. By adiabatically 
ramping through the resonant region, the quantum state 
will strongly overlap with the transporting state C, thus 
leading to an increase of the current. 

We compute the time-dependent current as j = 



1 



_ ^ ^ J p dt with p = {tp\p\tlj) which becomes the 

asymptotic ciirrent J in the limit t ^ oo. Fig. 2 displays 
the current evolution as we ramp with different rates. We 



find that the final current increases by an order of mag- 
nitude when the ramping velocity is lowered from the di- 
abatic into the adiabatic regime. If we ramp again, but 
backwards, with the same adiabatic rate, we will more or 
less return to the original state and momentum. However, 
if we ramp fast (diabatically) backwards, we will perform 
a Landau-Zener transition, and arrive at the initial value 
6 = —1.1, but now with a much larger current, than the 
one we started with! Fig. 2 shows the corresponding nu- 
merical results, as we ramp forward and backward with 
different ramping velocities. We ramp forward and back- 
ward between 6 = —1.1 and 6 = —0.9. When ramping 
forward we use a = 10~^ which corresponds to an adia- 
batic regime. During the backward ramping we use two 
values of a. Whereas in an adiabatic backward ramp the 
current becomes small again, for the diabatic ramp the 
current remains around its maximum value [Fig. 2]. 

The gap of an avoided crossing determines the maxi- 
mum adiabatic ramping rate. In order to increase that 
rate, the avoided crossing gap has to be increased. At the 
same time, for a Landau-Zener description to be applica- 
ble, the gap value should remain considerably smaller than 
the average level spacing. The gap value is determined by 
the strength of interaction between Floquet states. The 
interaction strength in turn is the larger, the more asym- 
metric the Floquet operator is. Motivated by this obser- 
vation as well as the work of Ref. [14], we further con- 
sider a bichromatic optical lattice potential to enhance 
the gap value. In particular, we consider the potential 
v{x) = cos(x) -|-ecos(2a;-|-0j,) as in Ref. [14]. We compute 
the quasienergy spectrum at ^ = —0.99 and ramp e. Fig. 
3 shows two avoided crossings of different states which lead 
to peaks with opposite current signs. The corresponding 
Husimi functions (see [11, 12] for details) reveal that such 
peaks appear due to resonances between a state in the 
chaotic layer and transporting states that carry opposite 
mean momenta. 

We use that structure and ramp adiabatically through 
this parameter window. As a result the initial state I with 
nonzero momentum transforms into a state with opposite 
momentum, and thus we reverse the current (a2 denotes 
the ramping rate for e). 

We first ramp e from —0.11 to —0.085 with state I as ini- 
tial input. After that, we compute the momentum of the 
obtained intermediate state as a function of time. Then 
we ramp e from —0.085 to —0.065, with the intermediate 
state as initial input. Then, we compute the evolution of 
the momentum of the final state. 

Fig. 4 shows the momenta of the intermediate and final 
states as a function of time. Both show good agreement 
with the momenta of the states II and III of Fig. 3, respec- 
tively. The oscillations indicate interference due to the 
coupling between the states. Husimi functions for states 
in adiabatic ramping show a good correspondence to the 
Floquet states II and III in Fig.3. We therefore manage 
to reverse and increase the current value significantly. 

If ramping 6 instead within the interval (—1.5,-0.5) 
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Fig. 3: Upper panel: (a) Quasienergy e vs amplitude e of the 
second harmonic in the potential v{x). The arrows indicate the 
ramping process starting in I, then passing through the inter- 
mediate state II, finally ending up in III. Inset: Quasienergy vs 
6 at the resonance point e = —0.097. (b) Current J vs e. The 
arrow indicates a narrow peak which appears due to an avoided 
crossing not shown in the graphic. Lower panel: Husimi func- 
tions of the states labeled by Roman numbers in upper panel 
(a). The intensity scales, in descending order, from yellow to 
black. We use the Husimi function defined in [12] . The param- 
eters are = —0.99 and Op = — 7r/2. The other parameters are 
the same as in Fig.l. 



Fig. 4: (a) Momenta of the intermediate and final states while 
ramping e as a function of time depicted by thick curves. The 
thickness of the curves is a result of periodic oscillations of 
the momenta with period T, which are not resolved by the 
eye in this time scale. Mean momentum for the states I, II 
and III in Fig. 3 are depicted by dashed lines. The dashed line 
for the state II appears superimposed to the momentum of 
the intermediate state. The ramping rate is Q2 ~ 2.5 x 10~^. 
(b) Husimi functions for intermediate and final states at t = 
50998 T. The other parameters are the same as in Fig. 3. 



(see inset of Fig. 3a), we find that, for an adiabatic ramp- 
ing, the transition rate should be below a = 4x 10^^. This 
value is much larger than the one for the symmetric case 
e = 0. Certainly, further increase of the gap values can 
be obtained by tuning other system parameters as well, 
thus further reducing the time needed for the adiabatic 
ramping. 

To summarize, we ramped control parameters through 
avoided crossings of Floquet states with different aver- 
age current characteristics. When diabatically ramping, 
the crossing is barely visible due to Landau-Zener transi- 
tions. For adiabatic ramping, the final quantum state will 
strongly change its transport properties. One can there- 



fore switch the average current from small to large values 
or vice versa. By first ramping adiabatically, and then 
ramping back diabatically, we can therefore reach a com- 
pletely different quantum state at one and the same value 
of the control parameter. Harvesting on several consecu- 
tive crossings, it is possible to achieve current reversals. 
In addition, by considering a bichromatic optical lattice 
potential, we show it is possible to increase the ramping 
speed while maintaining adiabaticity. 
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